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This paper discusses supersymmetry discovery potential of the upgraded 
D0 experiment at the Tevatron pp collider. Six final states with large trans- 
verse energy (momentum) leptons or photons (with or without large transverse 
momentum imbalances) are studied. These final states are expected to have 
small backgrounds and are thereby ideal for new physics searches. Implica- 
tions of the analyses of these final states on Gauge Mediated Supersymmetry 
Breaking models are discussed for integrated luminosities of 2 and 30 fb _1 . 
This study demonstrates that a large class of supersymmetry models can be 
probed at the upgraded Tevatron. 
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I. INTRODUCTION 



This paper summarizes the work done for the Tevatron Run II Higgs/Supersymmetry 
workshop |L| on the supersymmetry models with Gauge Mediation/Low Scale Supersymme- 
try Breaking (GMSB) ||. Six final states in which new physics might manifest itself are 
investigated using the parameters of the upgraded D0 detector ||. All of these final states 
are expected to have small physics and instrumentation backgrounds. Implications of the 
analyses of these final states in future Tevatron runs on the minimal (and not-so-minimal) 
GMSB models are discussed. Estimated discovery reaches in the supersymmetry parameter 
space are presented. 

II. OBJECT IDENTIFICATION 

Due to a large number of Monte Carlo (MC) events generated, no detector simulation 
is done for supersymmetry signals. All studies described in this paper except those ex- 
trapolated from Run I analyses are carried out at the particle level of the Isajet MC 
program [|J. A 2 TeV Tevatron center-of-mass energy is assumed throughout the studies. 
Leptons (£ = e, /i) and photons (7) are 'reconstructed' from the generated particle list by 
requiring them to have transverse energy (Ex) or momentum (p<r) greater than 5 GeV and 
to be within the pseudorapidity ranges: 

e: < 1.1 or 1.5 < \rj\ < 2.0; 

/z: \r]\ < 1.7; 

7: < 1.1 or 1.5 < \r]\ < 2.0. 

These fiducial ranges are dictated by the coverages of the electromagnetic calorimeter and 
the central tracker of the D0 detector. Furthermore, the leptons and photons must be 
isolated. Additional energy in a cone with a radius 1Z = J (A0) 2 + (A77) 2 = 0.5 in r\ — 
space around the lepton/photon is required to be less than 20% of its energy. 

Jets are reconstructed using a cone algorithm with a radius 1Z = 0.5 in rj — space 
and are required to have E T > 20 GeV and \if\ < 2.0. All particles except neutrinos, the 
lightest supersymmetric particles (lsp), and the identified leptons and photons are used in 
the jet reconstruction. The transverse momentum imbalance (-^t) is defined to be the total 
transverse energy of neutrinos and the LSPs. 

Energies or momenta of leptons, photons and jets of Monte Carlo events are taken from 
their particle level values without any detector effect. Smearing of energies or momenta 
of leptons, photons and jets according to their expected resolution typically changes signal 
efficiencies by less than 10% relatively and therefore has negligible effect on the study. 

The reconstruction efficiencies are assumed to be 90% for leptons and photons. For the 
purpose of background estimations, the probability (V) for a jet to be misidentified as a 
lepton (j — > £) or a photon (j — > 7) is assumed to be 10~ 4 . The probability for an electron 
to be misidentified as a photon (e — > 7) is assumed to be 4 x 10~ 4 , These probabilities are 
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a factor of three or more smaller than those obtained in Run I []. With a new magnetic 
central tracking system and a fine-segmented preshower detector, they should be achievable 
in Run II. 

In Run I, tagging of b-jets was limited to the use of soft muons in D0. Secondary vertex 
tagging of b-jets will be a powerful addition in Run II. For the studies described below, a 
tagging efficiency of 60% is assumed for those b-jets with Ej> > 20 GeV and \rj\ < 2.0. The 
probability V(j — > b) for a light-quark or gluon jet to be tagged as a b-jet is assumed to be 
10~ 3 . These numbers are optimistic extrapolations of what CDF achieved in Run I. 

Heavy stable charged particles can be identified using their expected large ioniza- 
tion energy losses (dE/dx) in the silicon detector, fiber tracker, preshower detectors and 
calorimeter. Based on Ref. f|, a generic dE/dx cut is introduced with an efficiency of 68% 
for heavy stable charged particles and a rejection factor of 10 for the minimum ionization 
particles (MIP). Note that the efficiency for identifying at least one such particle in events 
with two heavy stable charged particles is 90%. 

With the addition of preshower detectors, D0 will be able to reconstruct the distance 
of the closest approach (dca) of a photon with a resolution ~ 1.5 cm ||. Here the dca 
is defined as the distance between the primary event vertex and the reconstructed photon 
direction. Thereby it will enable us to identify photons produced at secondary vertices. In 
the following, a photon is called displaced if its dca is greater than 5.0 cm and is denoted 
by 7'. We further assume that the probability for a photon produced at the primary vertex 
to have the measured dca> 5 cm is P(7 — > V) = 2 x 10- 3 (about 3a). 

All final states studied have large Et (pt) leptons/photons with or without large Et- 
Triggering on these events are not expected to have any problem. Nevertheless, we assume 
a 90% trigger efficiency for all the final states. 



*The typical numbers determined in Run I are V(j — > e) = 5 x 10 4 , V(j — > 7) = 7 x 10 4 , and 
P(e^7) = 4x 10- 3 . 



III. FINAL STATES 



Signatures for supersymmetry vary dramatically from one model to another. They can 
also be very different for different regions of the parameter space in a model. Furthermore, 
these signatures are generally not unique to supersymmetry. In fact, some of the signatures 
are also expected from other theories beyond the Standard Model. Instead of chasing af- 
ter theoretical models (all of which except perhaps one are wrong anyway) a set of final 
states which are somewhat generic to many new physics models, including supersymmetric 
models, is identified. All of these final states are characterized by high Et{pt) isolated 
leptons/photons with or without large missing transverse momentum and are thus expected 
to have small physics and instrumental backgrounds. In the following, we discuss selection 
criteria and estimate observable background cross sections for six such final states: 



77-^T di-photon events with large Ifo 

jbjlftr single-photon events with b-jets and large I^t 

j'jjtfx single displaced photon events with jets and large I£t 

it + dE jdx high px di-lepton events with large ionization energy 

llljlpT tri-lepton events with jets and 

i^i^jjl^T like-sign di-lepton events with jets and I$t 



A. Final State 

The D0 Collaboration reported a search for di-photon events with large Ifx 
events) motivated by supersymmetric models with a light gravitino (G) as the lsp from a 
data sample of an integrated luminosity (£) of 106.3 ± 5.6 pb _1 in Run I. The 77-^r events 
were selected by requiring two identified photons, one with Ej, > 20 GeV and the other with 
Ej* > 12 GeV, each within pseudorapidity |r] 7 | < 1.1 or 1.5 < |?7 7 | < 2.0, and a ffiT greater 
than 25 GeV. Two events satisfied all requirements. 

The principal backgrounds were multijet, direct photon, W + 7, W + jets, Z — > ee, and 
Z — > tt — > ee events from Standard Model processes with misidentified photons and/or 
mismeasured Ifa. The numbers of estimated background events were 2.1 ± 0.9 from If T 
mismeasurement (QCD) and 0.2 ± 0.1 from misidentified photons (fakes). This led to an 
observed background cross section of 20 fb from QCD and of 2 fb from fakes in Run I. The 
I$t distributions before the I$t cut for both candidates and background events are shown in 
Fig. [I]. Note that events with large I$t are rare. 

Since the backgrounds are dominated by the I$t mismeasurement, they can be signifi- 
cantly reduced by raising the I^t cut. Therefore, the following selection criteria are used for 
the Run II studies: 

1) At least two photons with Ej, > 20 GeV; 

2) E T > 50 GeV. 
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The backgrounds with this set of selection criteria are expected to be significantly reduced 
by the increased cutoffs on Ifi? and photon and by the improved photon identification. 
The total observable background cross section in Run II is estimated to be cx 6 = 0.4(QCD) + 
0.2(fakes) = 0.6 fb assuming reduction factors of 5 from the raised IpT cutoff, 4 from the 
improved V(j — * 7), 3 from the higher photon Et requirement, and 10 from the decreased 
e — ► 7 fake probability. 
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FIG. 1. The If? distributions of the 77 and background samples for Run I. The number of events 
with ]$t< 20 GeV in the background sample is normalized to that in the 77 sample. Note that 
there is a 14 GeV I$t requirement in the trigger. The Ifix values plotted here are calculated off-line 
and therefore may differ from their values at the trigger level. 




B. -/bj^T Final State 

D0 Collaboration carried out a search M for single-photon events with at least two jets 
and large I£t {iJJ^t events) in Run I. The jjjtfr events were selected by requiring at least 
one identified photon with Ej, > 20 GeV and within pseudorapidity ranges |^ 7 | < 1.1 or 
1.5 < |?7 7 | < 2.0, two or more jets having E J T > 20 GeV and \r] j \ < 2.0, and If T > 25 GeV. 
A total of 318 events were selected from a data sample corresponding to an integrated 
luminosity of 99.4 ± 5.4 pb _1 . 

The principal backgrounds were found to be QCD direct photon and multijet events, 
where there was mismeasured I^t and a real or fake photon. The number of events from 
this source was estimated to be 315 ± 30. Other backgrounds such as those from W with 
electrons misidentified as photons were found to be small, contributing 5 ± 1 events. This 
led to an observed background cross section of 3,200 fb from the IpT mismeasurement and of 
50 fb from the fakes. The Ifix distribution before the Ifr> 25 GeV cut is shown in Fig. ^. As 
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shown in the figure, the backgrounds can be significantly reduced by raising the requirement 
on If}?- 

Events with a high Et photon, b-jets and large Ifi? are expected in some new physics 
models. These events, referred to as ^bjlpT-, are in many ways similar to the ^fjj^r events 
and thereby can be selected similarly: 

1) At least one photon with E? > 20 GeV; 

2) At least two jets with E 3 T > 20 GeV; 

3) At least one jet is tagged as a b-quark jet with E\ > 20 GeV; 

4) E T > 50 GeV; 

5) No leptons with E% > 20 GeV. 

The backgrounds from the QCD multijet events with real or misidentified photons and 
from the W events with electrons faking photons are estimated to be 0.63 fb, assuming 
background reduction factors of 5 from the raised Ifix requirement, 2 from the improved 
photon identification and using the assumed value of V(j — > b). The dominant background 
sources are expected to be ■jbb and ■yti events. These background sources cannot be reduced 
by the tagging of b-jets. However, the ■jbb contribution is expected to be small due to 
the large Ifi? requirement. Monte Carlo studies show that it is negligible. The jti (with 
tt — > W + W~bb) contribution is reduced by the requirements 4) and 5) and is estimated 
using the cross section of Ref. 0. A total of 0.9 fb observable background cross section is 
assumed. 

C. ijitfr Final State 

Photons produced at secondary vertices are predicted in a class of new physics models. 
These photons will appear to have large values of dca. Though dramatic, they alone are 
unlikely to be sufficient to reduce backgrounds from cosmic rays or from mismeasurement. 
We therefore select events with displaced photons accompanied by jets and large tfr'- 

1) At least one displaced photon with EP^ > 20 GeV; 

2) At least two jets with E 3 T > 20 GeV; 

3) E T > 50 GeV. 

These are called j'jj^T events. The dominant backgrounds are the same as those for the 
jjjl^T events, with a vertex-pointing photon being misidentified as a displaced photon. 
Using 7 7 (7 — > 7'), the observable background cross section from QCD and W events is 
estimated to be 0.6 fb. 
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FIG. 2. The I$t distributions of the jjj and background events. The number of events in the 
background is normalized to the jjj sample for I/t < 20 GeV, the region left of the dot-dashed 
line. 



D. High p T U + dE/dx Final State 



One possible new physics signature is the presence of heavy stable charged particles. 
These particles, if produced, will manifest themself in the detector as slowly moving muons 
with large ionization energy losses. Though D0 had several di-lepton analyses in Run I, 
none of these can be extrapolated to Run II, thanks to the replacement of the central 
tracker. Based on the expected signatures of several supersymmetric models with heavy 
stable charged particles discussed below, we select high pt di-lepton events (££ + dE/dx) 
with large dE/dx loss using the following requirements: 

1) At least two leptons with p l T > 50 GeV; 

2) M u > 150 GeV; 

3) At least one lepton passing the dE/dx requirement. 

The di-lepton mass requirement is intended to reduce Drell-Yan backgrounds. The principal 
backgrounds are: QCD dijet events with jets misidentified as leptons, tt, and Drell-Yan 
events. Using V(j — > £) — 10~ 4 and the assumed rejection factor of the dE/dx cut for 
the MIP particles, the observable background cross sections are estimated to be 0.1 fb from 
QCD dijet, 0.2 fb from tt events, and 0.2 fb from Drell-Yan processes. The QCD dijet cross 
section for px > 50 GeV is assumed to be 1 pb in the estimation. The total observable cross 
section is therefore 0.5 fb for the above selection. 
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E. mjtf T Final State 



D0 searched for gaugino pair production using the tri-lepton signature |TD[ in Run I. The 



lepton p T cut was typically 15 GeV for the leading lepton and 5 GeV for the non-leading 
leptons. The analysis also had a small I^t requirement. The observable background cross 
section was estimated to be around 13 fb. Most of these backgrounds are due to Drell-Yan 
processes. We select the £££jl$T events using the following criteria: 

1) Vt > 15 GeV > Pt > 5 GeV > Pt > 5 GeV ; 

2) E T > 20 GeV; 

3) At least one jet with E 3 T > 20 GeV. 

The Drell-Yan production, a major background source for the Run I analysis, is significantly 
reduced by the new jet requirement. The total observable background cross section is esti- 
mated to be 0.3 fb assuming background reduction factors of 10 from the jet requirement, 2 
from the improved particle identification, 2 from the higher If? cut. 



F. ^l^jj^T Final State 

Like-sign di-lepton events are expected from processes such as gluino pair production. 
They are also expected from processes with three or more leptons in the final states, but only 
two are identified. This final state is expected to have small backgrounds. Again without a 
magnetic tracker, D0 had no analysis of this nature in Run I. Based on Monte Carlo studies 
for several supersymmetric models, we select i^i^jjl^T events using the following criteria: 

1) Two like-sign leptons with > 15 GeV; 

2) At least two jets with E T > 20 GeV; 

3) E T > 25 GeV. 

Events with three or more identified leptons are removed to make the sample orthogonal 
to the iiijl^T sample. Since leptons are relatively soft in pj- for the new physics model 
we investigated using this selection, the effect of charge confusion due to a limited tracking 
resolution is thus neglected in this study. The major backgrounds are: W + jets events 
with one of the jets misidentified as a lepton, ti events with energetic leptons from b-quark 
decays, and Drell-Yan (WZ, ZZ) events. The W + jets background is estimated using the 
number of W + 3j events observed in Run I folded with V(j — > £) to be 0.2 fb. The ti and 
Drell-Yan backgrounds are estimated using Monte Carlo to be 0.1 and 0.1 fb respectively. 
Adding the three background sources together yields a total observable background cross 
section of 0.4 fb. 
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IV. CONSTRAINTS ON GMSB MODELS 



The supersymmetric models with gauge mediated supersymmetry breaking are charac- 
terized by a supersymmetry breaking scale A as low as 100 TeV and a light gravitino which is 
naturally the lightest supersymmetric particle. In these models, supersymmetry is assumed 
to be broken in a hidden sector and the symmetry breaking is transmitted to the visible sec- 
tor of Standard Model particles and their superpartners through the Standard Model gauge 
interactions. The minimum gauge mediated supersymmetry breaking model is described by 
five parameters: 



A: supersymmetry breaking scale 

M m : messenger sector scale 

N: number of messengers 

tan (3: ratio of the v.e.v of the two higgs doublets 

sign(/i): sign of the higgs mass parameter 



The phenomenology is largely determined by the next-to-lightest supersymmetric particle 
which in turn depends on the values of the above five parameters. For a review of GMSB 
models, see Ref. |§. In the following, we discuss expected sensitivities with integrated 
luminosities of 2, 30 fb _1 for four different model lines defined by the working group. Each 
model line has a different nlsp. All theoretical expectations and signal efficiencies are 
obtained from the Isajet MC program. A minimum px of 50 GeV of the hard scattering 
is applied for all signal processes. We define the significance (N s /SN b ) as the ratio between 
the number (N s ) of expected signal events and the error (8N b ) on the estimated number 
of background events. Here a 20% systematic uncertainty is assumed for all estimated 
observable background cross sections. Therefore, 

5N b = ^/C-a b + (0.2 • £ ■ a b f 

We characterize the sensitivity using the minimum signal cross section <Jdis for a 5 standard 
deviation (5<r) discovery: 

N s _ £ ■ o d is ■ e _ _ 
5N b ~ 6N b 

where e is the efficiency for the signal. The minimum observable signal cross section <J obs 
defined as er^is ' e f° r the discovery is therefore independent of signal processes. The a obs as a 
function of £ for several different values of a b are shown in Fig. |3|. It decreases dramatically 
as £ increases for small £ values and flattens out for large £ values. Clearly, the sensitivity 
can be improved for large £ values by tightening the cuts to reduce backgrounds further. 

In the following, we express the 5<r discovery cross sections as functions of the super- 
symmetry breaking scale A and the lighter chargino (xt) mass for the four different model 
lines. 
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FIG. 3. The minimum observable signal cross section a \, s for a 5a discovery as a function of 
integrated luminosity for four different values of background cross sections. 



A. Model Line 1: Xi as the NLSP 



Within the framework of the minimal GMSB models, Xi is the NLSP for most of the 
parameter space. If the Xi has a non-zero photino component, it is unstable and decays to 
a photon plus a gravitino (xi ~^ 7^0 with a branching ratio of nearly 100%. Depending 
on its lifetime, pair production of supersymmetric particles will result in 7-^t, and 

Ifr+X events. For the purpose of this study, we consider a class of models with the following 
parameters fixed: 

N=l, ^ = 2, tan/? = 2.5, fj, > 

while A is allowed to vary. For the range of A values of interest at the Tevatron, the super- 
symmetry production cross section is dominated by XiXi an d X1X2 production. Figure [| 
shows the schematic decay chains of Xi an d X® with their branching ratios for A = 100 TeV. 
In the following, scenarios with prompt and delayed Xi decays are discussed. If the Xi is 
quasi-stable, i.e. with a long lifetime, the signature will be identical to that of the super- 
symmetric models with gravity mediation. 



1. Prompt X\ —> jG Decay 

If the xi decays in the vicinity of the production vertex, events are expected. 

The distributions of photon Et and event I^t for A = 80, 140 TeV are shown in Fig. |5|. 
These events typically have high Et photons together with large transverse momentum 
imbalances, and therefore can be selected using the 77-^T criteria discussed in Section HI A . 
Table | shows the detection efficiencies, significances along with the total theoretical cross 
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FIG. 4. Decay schematics of Xi an d x\ f° r A = 100 TeV for the model line with a x? as the 
nlsp. Percentage branching ratios for main decay modes are shown in parentheses. 



sections of super symmetry, chargino and neutralino masses for different values of A. Figure |6| 
compares the 5cr discovery cross sections a^ s with the theoretical cross sections expected 
from supersymmetry for two different values of £ as functions of the lighter chargino mass 
m~± (and the supersymmetry breaking scale A). The lighter chargino with mass up to 290, 
340 GeV can be discovered for C=2, 30 fb" 1 respectively. 



A (TeV) 


60 


80 


100 


120 


140 


160 


oth (fb) 


464 


105 


27 


7.7 


2.2 


0.7 


m-± (GeV) 


138 


194 


249 


304 


357 


410 


m^o (GeV) 


75 


104 


132 


160 


188 


216 




16.1 


24.3 


28.2 


30.1 


30.6 


30.2 


N s /5N b (2 fb" 1 ) 


136 


46 
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N s /5N b (30 fb" 1 ) 
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41 


12 


3.6 


1.1 



TABLE I. The supersymmetry cross section (a t h), Xi an d Xi masses, detection efficiency of 
the selection, and significances for different values of A for the models with a short-lived Xi 

as the NLSP (model line 1). The relative statistical error on the efficiency is typically 2%. The 
observable background cross section is assumed to be 0.6 fb with a 20% systematic uncertainty. 



2. Delayed xi —* lG Decay 



If the Xi has a significant lifetime, the photon from its decay may not point back to 
the primary vertex. If the decay occurs inside the tracking volume of the D0 detector, the 
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FIG. 6. The 5a discovery cross section curves as functions of mass of the lighter chargino (and 
the supersymmetry breaking scale A) along with the theoretical cross sections for the model line 
1. The nlsp Xi is assumed to be short-lived. Two curves corresponding to integrated luminosities 
of 2 and 30 fb -1 are shown. 
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FIG. 7. (a) Average decay distance and DCA for different values of proper decay distance cr 
of the nlsp Xi- (b) The probability that a photon decays within the tracking volume and the 
probability that such a photon has its DCA> 5 cm as functions of cr. Note that A = 100 TeV 
(m~± = 249 GeV) for both plots. 



photon is expected to traverse standard electromagnetic detectors (the preshower detectors 
and the electromagnetic calorimeter). It, therefore, can be identified. However, if the decay 
occurs outside the tracking detector, the photon identification is problematic. For this study, 
we assume that the photon is identifiable if it is produced inside a cylinder denned by the 
D0 tracking volume (r < 50 cm and \z\ < 120 cm) and is lost if it is produced outside the 
cylinder. Figure [F|(a) shows the average decay distance and distance of closest approach of 
the X? as functions of its proper decay length (cr) for A = 100 TeV. Due to its heavy mass, 
the Lorentz boost for the Xi is typically small (7 ~ 1.5). The probabilities that a photon 
is identifiable and that an identifiable photon has dca> 5 cm as functions of the Xi proper 
decay distance cr are shown in Fig. 0(b), again for A = 100 TeV. Distributions for other A 
values are similar. Figure |] shows jet multiplicity and Et distributions for A = 80, 140 TeV. 
Most of these events have large E T jets and thus can be selected using the j'jjtfr selection 
criteria discussed above. The detection efficiencies and the expected significances of the 
Yjjl^T selection for cr = 50 cm are tabulated in Table |1| as an example. The estimated 
5cr discovery reaches in A and chargino mass for different values of cr are shown in Fig. |9| 
along with those expected from the analysis. As expected, the 77-^r analysis has a 

stronger dependence on cr than the j'jj^T analysis. 



B. Model Line 2: f x as the NLSP 



If fi (the lighter of the two mixed states of tr and tl) is lighter than Xii au super- 
symmetric particles will cascade into the f\ which in turn will decay to tG with a 100% 
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(m-± = 194, 357 GeV) for the models with Xi as the nlsp. All distributions are normalized 
to have unit area. 
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TABLE II. The detection efficiency of the j'jj^T selection, and the significances for different 
values of A for the models with Xi as the nlsp. The Xi proper decay distance cr is assumed to be 
50 cm. The relative statistical error on the efficiency is typically 2%. The observable background 
cross section is assumed to be 0.6 fb with a 20% systematic uncertainty. 
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FIG. 9. The 5a discovery reaches of the 7' ' jjI$T an d the 77-$r analyses in m-± and A as functions 
of the proper decay distance of the Xi NLSP for C=2, 30 fb _1 . 



branching ratio. This class of models is defined by following parameter values: 

N = 2, = 3, tan/3 = 15, /i > 

with varying A. Again, xf xf an d xf X2 dominate the production cross section for A < 
75 TeV. For A values above 75 TeV, fifi, e^eR, and /i_R/i_R productions become important. 
As an example, branching ratios of Xi an d X2 f° r A = 40 TeV are graphically displayed 
in Fig. [31]. In the following, two cases corresponding to short-lived and quasi-stable f]S are 
discussed. It should be noted that the two analyses discussed below are also sensitive to the 
case with a intermediate f\ lifetime. 



1. Prompt r\ — > tG Decay 

If the f\ is short-lived and decays in the vicinity of the production vertex (i.e. with a 
decay distance 7cr < 10 cm), anomalous r production is expected from supersymmetry. 
Together with the W*/Z* productions from the cascade decays of primary supersymmetric 
particles, these events will give rise to £££j^T and £ ± £ ± jj^'t final states. The lepton pr 



distributions of the £££j^T and t t jjlfix events are shown in Fig. |TT]. Since most leptons 



are produced in r decays, their pxs are relatively soft. Table |T| shows the efficiencies of 
the iiijl^T an d t I jjlfc selection criteria for these events along with the theoretical cross 
sections, Xi an d ?\ masses. Note that the £££jI$T and £ ± £ ± jj-^T criteria are orthogonal. 
The efficiencies are relatively small largely due to the small branching ratio of the events to 
tri-leptons. We note that the total efficiencies shown in the table are somewhat conservative. 
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FIG. 10. Decay schematics of X\ an d X2 f° r A = 40 TeV for the model line with a ti as the 
nlsp. Percentage branching ratios for main decay modes are shown in parentheses. 



t events to the i^l^jj^T events 



They do not take into account the migration of the 
due to inefficiency in the lepton identification. The 5a discovery curves are shown in Fig. |I2 . 
The lighter chargino with mass up to 160 and 230 GeV can be discovered for C=2, 30 fb -1 . 

The conventional wisdom is that this analysis should benefit from a r identification. 
However, we doubt that it will have a dramatic impact on the reach in the supersymmetry 
parameter space. Though a r identification could improve the efficiency for the signal, it 
will undoubtably come with large backgrounds. Nevertheless, a r identification is essential 
to narrow down theoretical models if an excess is observed in the tri-lepton final state. 



2. Quasi-stable f\ 

If the f\ has a long lifetime (quasi-stable) and decays outside the detector (7CT greater 
than ~ 3 m), it will appear in the detector like a muon with the exception of a large 
ionization energy loss. The signature is, therefore, two high p? 'muons' with large dE/dx 
values. These events can be selected using the criteria described in Section [III D| . The 



expected pt distributions of the f\ for two different values of A are shown in Fig. [T3"|(a). 
The cut of pt > 50 GeV of the ££ + dE/dx selection is efficient for the signal while it is 
expected to reduce backgrounds significantly. The typical invariant mass of the two 'muons' 
(assuming massless) is very large as shown in Fig. 0(b). A Mu > 150 GeV requirement 
does little harm to the signals. Due to its large mass, the f\ is expected to move slowly. 
However since most of the fi's are produced in the decays of massive xfs and X2 S > t ne average 
speed /?(= v/c) is relatively large. It is around 0.7 for the A values studied. Note that the (3 



distribution is very similar to that shown in Fig. |T8|(b) for the models with I as the Co-nlsp 
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FIG. 11. Lepton px distributions for (a) the iiij^x events and (b) the l^i^jjlfir events of the 
models with a short-lived f\ for m~± = 182 GeV (A = 40 TeV). All distributions are normalized 

Xi 

to unit area. 
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TABLE III. The theoretical cross section, Xi an d f\ masses, detection efficiency of I I jjtfr 
and tttjI$T selections, and significances for different values of A for the models with a shortdived 
fi as the nlsp. The relative statistical error on the efficiency is typically 25%. The combined 
^^jjtfr and UtjI$T background cross section is assumed to be 0.7 fb with a 20% systematic 
uncertainty. 
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FIG. 12. The 5a discovery cross section curves as functions of mass of the lighter chargino and 
the supersymmetry breaking scale A for the model line 2 along with the theoretical cross sections. 
The 5ct curves are shown for both short-lived and quasi-stable nlsp's and for integrated luminosities 
of 2, 30 fb -1 . 
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FIG. 13. The lepton (a) and di-lepton mass Mm (b) distributions for the models with a 
quasi-stable f\ as the nlsp for A = 40, 80 TeV (m-± = 182, 394 GeV). The vertical dashed lines 
indicate the cutoffs. All distributions are normalized to unit area. 
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FIG. 14. The I$t distributions for the models with a quasi-stable f% as the nlsp for 

A = 40, 80 TeV (m-± = 182, 394 GeV). The decays xt -> X\W ± -> f 1 rW ± and X \ e e, M/* 

1 

contribute to events with small $r while the decays xt ~* an d X% — * ttj are the source for 
events with large 



Nevertheless, the not-so-slow moving fi's are expected to deposit large ionization energies in 
the detector, differentiating them from other high p? MIP particles. Since the backgrounds 
for the requirements px > 50 GeV and M« > 150 GeV are already small, it pays to have a 
dE/dx requirement with a relatively high efficiency for the signal and a resonable rejection 
for the MIP particles. The Ifix distribution of these events as shown in Fig. [H] shows two 
distinct regions: small and large Ifix- The decays xt ~^ XiW — > f\rW and X 2 ~^ f^fi 
contribute to events with small Ifa>. The decays xt ~^ an d X2 ~ * are responsible 
for events with large IpT- The detection efficiencies and the expected significances of the 
II + dE/dx selection for different values of A are tabulated in Table [IV|. The high efficiency 
is largely due to the high momentum expected for the quasi-stable f\. The 5a discovery 
curves are shown in Fig. 12 for two values of C. The lighter chargino with mass up to 340, 
410 GeV and the fi with mass up to 160, 200 GeV can be discovered for the two integrated 
luminosities respectively. 



C. Model Line 3: t as the Co-NLSP 



For some regions of the parameter space, three light sleptons (e^, flu, and fi) are es- 
sentially degenerate in mass and they can be lighter than Xv As a result, the sleptons 
(£ = fi, Sr, jlji) effectively share the role of the nlsp. Quantitative studies of this type of 
models are done for following GMSB parameter values: 

M m 

N = 3, — ^ = 3, tan (3 = 3, /1 > 
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TABLE IV. The theoretical cross section, \\ an d f\ masses, detection efficiency of the U+dE / dx 
selection, and significances for different values of A for the models with a quasi-stable f\ nlsp. The 
relative statistical error on the efficiency is typically 1%. The background cross section is assumed 
to be 0.5 fb with an uncertainty of 20%. 



with again varying A values. For small values of A, XiXi an d X1X2 dominate the production 
cross section. As shown in Fig. JSj, X1X1 an d X1X2 production will yield events with 
multileptons in the final state. The slepton pair production surpasses chargino-neutralino 
production if A > 50 TeV. The total supersymmetry cross section for several different A 
values can be found in Table [Vj The lifetime of I determines the event topology. In the 
following, we discuss the cases with short-lived and quasi-stable is. Again, the analyses 
should also be sensitive to the I nlsp with a intermediate lifetime. 
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FIG. 15. Decay schematics of X\ an d X2 f° r A = 40 TeV for the model line with scalar leptons 
as the Co-nlsp. Percentage branching ratios for main decay modes are shown in parentheses. 



20 



1. Prompt I — > iG Decay 

If the decay i — > £G is prompt (7cr < 10 cm), lllf^ events are expected from supersym- 
metry. Unfortunately, this final state is swamped by backgrounds from the Standard Model 
processes such as ti, WW, WZ and ZZ productions as well as from W + jets production 
with one of the jets misidentified as a lepton. However we note that these events typically 
have multiple leptons in the final state and most of them are in the central pseudorapidity 
region with good lepton identification. Apart from those from i decays, leptons are also 
expected from W*'s and Z*'s produced in the cascade decays of xt an d X2 °f supersym- 
metry originated events. Therefore, they can be selected using the iiijl^T criteria. The 
Pt distributions of the leading lepton and the third lepton of these events are shown in 
Fig. [TB|(a) . Since most of the leading leptons are produced in the direct decays of heavy 
£'s, its pt spectrum is relatively hard as shown in the figure. The detection efficiencies and 
the expected significances are summarized in Table [V]. The reduction in the relative cross 
section of the tri-lepton producing xtxt an d X1X2 processes is responsible for the decrease 
in efficiency as A increases. For A > 50 TeV, the it production cross section surpasses that 
of the xtxt an d xtx.2- With the i — > iG decay, ii events will result in a high p T iHfi T final 
state. We note that the improvement by adding the i I jjlpT selection is minimal in this 



case. The 5cr discovery curves are compared with the theoretical cross sections in Fig. [L7. 
With integrated luminosities of 2, and 30 fb" 1 , the lighter chargino with mass up to 310 and 
360 GeV can be discovered respectively. 
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FIG. 16. The pr distributions of (a) the leading lepton and (b) the third lepton for the models 
with a short-lived I as the Co-nlsp for m-± = 279, 437 GeV (A = 40, 60 TeV). Note that the 
Pt requirement is 15 GeV for the leading lepton and 5 GeV for the non-leading leptons. All 
distributions are normalized to unit area. 
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TABLE V. The theoretical cross section, \\ an d I masses, detection efficiency of the lltjlpT 
selection criteria, and significances for different values of A for the models with short-lived Ps as 
the Co-nlsp. The efficiencies typically have a relative statistical uncertainty of 4%. The observable 
background cross section is assumed to be 0.3 fb with a 20% systematic uncertainty. 
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FIG. 17. The 5a discovery cross section curves as functions of mass of the lighter chargino and 
the supersymmetry breaking scale A for the model line 3 along with the theoretical cross sections. 
The 5a curves are shown for both short-lived and quasi-stable t Co-nlsp's and for integrated 
luminosities of 2, 30 fb -1 . 
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P T (GeV) p 
FIG. 18. The lepton px (a) and the NLSP speed (3 (b) distributions expected for models with a 
quasi-stable £ as the Co-nlsp for m~± = 279, 437 GeV (A = 40, 60 TeV), where f3 is measured in 
the unit of the speed of light c. The vertical dashed line indicates the cutoff. All distributions are 
normalized to unit area. 



2. Quasi-stable I 

If the £ has a long lifetime, it can decay outside the detector (jct > 3 m). In this case, 
the £ will appear in the detector like a 'muon' except that the ionization energy loss will 
be large. This signature is identical to that of a quasi-stable fi discussed above. Therefore, 
the signal events can be identified using the ££ + dE/dx selection. The expected pr and 



P distributions of the £ for A = 40, 60 TeV are shown in Fig. [18|. The £s typically have 



very large pt and are mostly central. For example, about 90% of the £s are in central 



pseudorapidity region with the tracking coverage for the case of A = 70 TeV. Table VI 



shows the detection efficiencies and the expected significances for different A values. The 5<r 



discovery curves are shown in Fig. 17. The lighter chargino mass discovery reach is about 



390 GeV for C=2 fb" 1 and 480 GeV for £=30 ftr 1 . 



D. Model Line 4: h as the NLSP 

For most of the parameter space, Xi wu l predominantly decay to 7G if it is the NLSP. 
However if Xi is higgsino-like (h), Xi ~~ * an d Xi ~^ hG decays could have significant 
branching ratios in some regions of the parameter space of non-minimal GMSB models. For 
the Run II studies, these models are defined to have fixed values of 

N = 2, -j^ = 3, tan/5 = 3, .A = -0.75 
A M 1 

with A varying. Here Mi is the mass parameter associated with the U(l)y symmetry. In 
these models, the lightest neutral higgs boson h has a mass around 104 GeV. Pair production 
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TABLE VI. The theoretical cross section, Xi an d £ masses, detection efficiency of the ££+dE/dx 
selection, and significances for different values of A for the models with quasi-stable ts as the 
Co-NLSP. The relative statistical error on the efficiency is typically 1%. The background cross 
section is assumed to be 0.5 fb with a systematic uncertainty of 20%. 



of supersymmetric particles may result in XiXi ~ * hG — ► ^bbGG final state which would 
appear as ybblpT events in the detector assuming prompt Xi decays. These events are 
characterized by high photons and large IpT as shown in the Fig. ^ for A = 80, 110 TeV. 
They can be selected using the jbj^x selection criteria discussed in Section |III B| . The 
detection efficiencies and N s /5Nb significances are shown in Table |V11| for different values 
of A. Most of the events selected are due to the 7/1 production with h — > bb. However, 
a non-negligible fraction of the events is actually due to the 'jZ + X with Z — > bb. The 
discovery reach in A and m-± is shown in Fig. |20] for £=2, 30 fb -1 . 




(GeV) E T (GeV) 

FIG. 19. The photon Et (a) and event J^t (b) distributions for models with a short-lived 
higgsino- liked xi a s the NLSP with m-± = 167, 229 GeV (A = 80, 110 TeV). All distributions are 
normalized to have unit area. 

If an excess is seen in the ^bjlfx final state, it will be of great interest to reconstruct the 
di-jet invariant mass. Figure [21] shows the invariant mass distribution of the two leading jets 
for A = 80 TeV and C=2 fb -1 . A mass peak around = 104 GeV is clearly identifiable. 
The asymmetry in the mass distribution is partly due to •yZ + X contribution and partly 
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FIG. 20. The 5a discovery cross section curves as functions of the super symmetry breaking scale 
A and the lighter chargino mass for the model line 4 along with the theoretical cross sections. The 
h is assumed to be short-lived. 
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TABLE VII. The theoretical cross section, \\ an d h masses, Xi decay branching ratios, detection 
efficiency of the jbjl^T selection, and significances for different values of A for the models with a 
short-lived higgsino as the NLSP. The relative statistical error on the efficiency is typically 4%. 
The background cross section is assumed to be 0.9 fb with a 20% systematic error. 
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due to the effect of gluon radiation and the energy outside the jet cone. 
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FIG. 21. The invariant mass distribution of the two leading jets for m-± 
(A = 80 TeV) and C=2 fb" 1 
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Note that 19 signal and less than 2 background events are ex- 



V. SUMMARY 

In this paper, observable background cross sections for the six final states in which new 
physics might manifest itself are estimated. All the final states studied are found to have 
small backgrounds. Implications of the analyses of these final states for future Tevatron 
runs are discussed in the framework of Gauge Mediated Supersymmetry Breaking models. 
Potential discovery reaches in supersymmetry parameter space for integrated luminosities of 
2 and 30 fb _1 are examined for models with different nlsp. Though the selection criteria are 
not optimized for the models discussed and not all final states are investigated, the study 
does show that the upgraded D0 experiment at the improved Tevatron collider has great 
potentials for discovery. 
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See the workshop web page: |http : //f rrth37 . f nal . gov/susy . html 
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